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Abstract. 

j^ . We present an efficient algorithm for computing the partition function of the q- 

jy>2 ' colouring problem (chromatic polynomial) on regular two-dimensional lattice strips. 

Our construction involves writing the transfer matrix as a product of sparse matrices, 
each of dimension ~ 3™, where m is the number of lattice spacings across the strip. 
As a specific application, we obtain the large-g series of the bulk, surface and corner 
'"^ . free energies of the chromatic polynomial. This extends the existing series for the 

Ch I square lattice by 32 terms, to order q~'^^. On the triangular lattice, we verify Baxter's 

analytical expression for the bulk free energy (to order q~^^), and we are able to 
conjecture exact product formulae for the surface and corner free energies. 
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S ' 1- Introduction 

iX ' The g-colouring problem consists in assigning to each of the vertices of a graph G any 

O I one of q different colours, in such a way that adjacent vertices carry different colours. 

The number of possible colourings (possibly zero) is known as the chromatic polynomial 
Pcil) of the graph. Although this definition supposes that g is a positive integer, it is 
K> ! not hard to show [H |2] that Pg{<i) is indeed a polynomial in g, and so the number of 

colours can be regarded as a formal variable. 

The history of the g-colouring problem is long and interesting, and we refer the 
reader to [3] for an extensive list of references. The case where G is a planar graph has 
attracted particular interest. A well-known result is then the four-colouring theorem |1] 
which states that -Pg(4) > for any planar G. Other results exploit the extension of 
Pg{(i) to non-integer values of q. One interesting question is whether there exists, in 
the planar case, some gc so that Pcil) > for all q G [qc, oo). This statement has been 
established as a theorem [5] for qc = 5, and the so-called Birkhoff-Lewis conjecture [5] 
that this extends to gc = 4 is widely believed to be true. 

Apart from these general results, the case where G is a regular planar lattice — 
typically square or triangular — has also been intensively studied. In particular, a long 
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series of papers [SI El El El 13 CHI [HI 112] is devoted to the study of the location and 
properties of chromatic zeroes, Pg{<1) = 0, in the complex g-plane for a variety of 
boundary conditions. Some of the mechanisms responsible for the generation of real 
chromatic zeroes in the region q G [0, 4] have even been understood analytically [T3|[Ti]. 
One should also mention that there exists a family of planar graphs with real chromatic 
roots converging to g = 4 from below [T5] . 

An important technical ingredient in the studies of chromatic zeroes is the 
possibility to build up Pcil), for regular lattice strips of a given width m, by means 
of a transfer matrix T. In order to elucidate the amazingly intricate behaviour of the 
chromatic polynomial in the thermodynamical limit, it is clearly desirable to be able to 
access as large m as possible in these computations. 

In the present paper we show how to improve the constructions of [Ml El El El E] by 
writing T as a product of sparse matrices, each of dimension ~ 3™. This is implemented 
in an efficient way, that enables us to compute exactly Pail) for rectangles as large as 
20 X 21. Using the finite-lattice method [IT] this results in series expansions — in powers 
of 1/q — of the bulk, surface and corner free energies of the chromatic polynomial. To 
give an idea of the improvement on existing results, the bulk free energy series on the 
square lattice was taken to order q~^^ in [18], and very recently extended to order q''^''^ 
in [11]. We here add a further 32 terms to this latter series, taking it to order g~^^. The 
surface and corner free energies are similarly obtained to orders g~^^ and g~^^. 

The case of the triangular lattice — which does not appear to have been studied 
previously in terms of series expansions of the free energy — turns out to be particularly 
interesting. Baxter [191 E] has shown that for this lattice the chromatic polynomial 
possesses an integrable structure. This permitted him to derive an analytical expression 
for the bulk free energy, in the form of an infinite product formula, valid in the region 
q G (— oo, 0] U [4, oo). We obtain in this case the bulk, surface and corner free energies 
to order g~'^°. Our results validate Baxter's result for the bulk free energy, and permit 
us to conjecture also the exact infinite product formulae for the surface and corner free 
energies. 

2. Chromatic polynomial 

We first recall how the chromatic polynomial appears as a special case of the partition 
function of the g-state Potts model. Given a finite undirected graph G = (V, E) with 
vertex set V and edge set E, the Hamiltonian of the Potts model is given by 

H = -jJ2 ^^..^. ' (1) 

where J is the coupling constant, and a spin cfi G {1, 2, . . . , g} is defined at each vertex 
i & V. The Kronecker delta function S„.^^. equals 1 if cij = Cj, and otherwise. 
Setting f = e'^ — 1, the partition function can be rewritten as 

ZG{q,v)=J2 n [1+^'^'^..-.] (2) 

{a} (Ji)eE 
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= J2 g'^(^)t;l^l , (3) 

ACE 

where k{A) is the number of connected components in the subgraph Ga = {V,^), and 
l^l is the number of elements in the edge subset A (1 E. In the rewriting ([3]) — known as 
the Fortuin-Kasteleyn (FK) representation |2^ — q does no longer need to be an integer, 
and in particular can be used as an expansion parameter. For g — )■ oo, only the term 
with A = — whence k{A) = \V\ — contributes, and thus one can develop Zc^q^v) in 
powers of 1/g. The chromatic polynomial appears as the special case v = —1, 

PGiq) = ZGiq,-l), (4) 

as is evident from ([2]). 

3. Transfer matrix formalism 

Thanks to the FK representation ([3]), the partition function for a lattice strip of width 
m can be built using a transfer matrix T of dimension independent of q fJGl [3] . In this 
section we briefly review the necessary formalism, with the minor modification that we 
rewrite things in terms of q~^ rather than q, since this is more convenient for our series 
expansion purposes. We then introduce several improvements: 

(i) A factorisation of T into sparse matrices that requires only the use of non-nearest- 
neighbour and almost non-nearest-neighbour partitions. 

(ii) An efficient labelling of non-nearest-neighbour partitions by natural numbers. 

(iii) The use of modular arithmetic instead of arbitrary-precision integer arithmetic. 

The combination of these ingredients will yield an efficient algorithm that allows, in 
particular, to add further terms to the series expansions of [iHl fTT] . 

Let us start by noting that, in view of (|3]), T can be taken to act on states which 
are partitions of the vertex set {1,2,..., m} corresponding to a single horizontal layer 
of spins. The role of T is to propagate the system one lattice spacing upwards in the 
vertical direction. The states describe how the connected components of Ga = (V, A) 
partition the m vertices of a layer, where A is the intersection of A with the set of 
edges which are below the layer under consideration. For instance, when m = 3, the five 
possible states are the partitions 

(1)(2)(3), (12)(3), (13)(2), (1)(23), (123). (5) 

We shall see shortly how the information contained in these states permits us to account 
for the non- locality of the Boltzmann weight q''^^\ 

Our objective is to write the chromatic polynomial on a rectangular piece, of width 
m and height n, of the square or triangular lattice in the form 

Pciq) = (eo|TiastT"~>o) (6) 

where the basis state |eo) corresponds to the all-singleton partition (1)(2) ■ ■ ■ (m). Each 
transfer matrix T adds one layer to the lattice, consisting of m — 1 horizontal edges and 



Free energy series for the chromatic polynomial 4 

m vertical edges — and for the triangular lattice also m — 1 diagonal edges. A special 
matrix Ti^st will be required to deal with the last layer of rra — 1 horizontal edges. The 
notation (eo| means projection on the all-singleton state, and is not meant to imply the 
definition of a scalar product on the space of partitions (although such a scalar product 
can indeed be defined). 

Let us observe that it would be quite wasteful to construct the matrix T explicitly. 
Not only would this require memory resources of the order of (dim T)^ for its storage, 
but computing the product T\v), where \v) represents an arbitrary linear combination 
of all possible basis states, would take a time proportional to (dim T)^. If, on the other 
hand, we write T as a product of matrices each corresponding to the addition of a single 
edge, these matrices will turn out to be sparse (with only one or two non-zero entries 
per column), and so both time and memory requirements are linear in the dimension of 
these matrices. In order to proceed, we first study these single-edge matrices in some 
detail. 

3.1. Single-edge transfer matrices 

Since we are going to construct a series expansion in powers of 1/g, it will be useful to 
pull out a trivial factor in ([3]) and study instead 

ZG{q,v) = q-^''^Za{q,v) = ^ g^'(^)-|V'lt;l^t . (7) 

ACE 

The lowest-order term, corresponding to A = 0, is then 1. It is useful to think of the 
factor g*''(^)~l^l in the following way: whenever the addition of an edge to the set A joins 
up two connected component the result is a factor of q'^. 

The single-edge transfer matrix that adds a horizontal edge between vertices i and 
i + 1 can now be written 

Hi,i+i = 1 + f Jj,i+i , (8) 

where the first (resp. second) term represents the possibility of that edge being absent 
from (resp. present in) the subset A. The join operator Jj^j+i acts on a partition of 
{1,2,..., m} as the identity if i and i + 1 belong to the same block. If they belong to 
two different blocks, the action of Jj j+i is to join up the blocks and attribute a factor 

We now make a remark that will turn out to be very important in what follows. 
For the case of the chromatic polynomial (i.e., v = —1) Hj j+i is idempotent. It projects 
out states in which i and 2 + 1 belong to the same block. It turns out useful to introduce 
the normalised projection operator Pi,i+i which acts on a partition as the identity if i 
and i + 1 belong to different blocks, and as zero if i and i + l are in the same block. We 
have then 

Hi,i+iPi,i+i = Hj^j+i = (Hj^j+i) . (9) 

The transfer matrix corresponding to a vertical edge between vertices i and i' can 
be thought of as the following sequence of operations: 
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(i) Insert i' into the partition as a singleton. 

(ii) Apply Hi,i/. 

(iii) Eliminate i from the partition. 

(iv) Relabel i' as i. 

Alternatively, the vertical edge as can be represented as the operator 

V, = vq~' + D, , (10) 

where the first (resp. second) term represents the situation where that edge is present 
in (resp. absent from) the subset A. The detach operator Dj acts by transforming i 
into a singleton. The advantage of the rewriting is that the relabelling z' — )■ i is done 
automatically, but one should keep in mind that i does not represent the same vertex 
on the lattice before and after the action of Vj. 

Note that although the vertical-edge transfer matrix Vj is related to Hj^j/, it is not 
a projector in the chromatic case v = —1. 

We next show how to combine the observation that Hj j+i is a projector with the 
sparse matrix decomposition of T in order to achieve a very efficient transfer matrix 
formalism. Although this can be done for any regular lattice, we shall discuss only the 
cases of the square and triangular lattices. 

3.2. Square lattice 

One can use the operators Hj j+i and Vj to build up the lattice, by writing T as a suitable 
product of such operators. An obvious construction is to first lay down all the horizontal 
edges of a layer, and the add all vertical edges that connect to the next layer: 

nV^jfnH.m) (anyt;). (11) 

Since the next transfer matrix in the product (E]) will first lay down a row of horizontal 
edges, we may use ([9]) and write instead, for t> = — 1, 

n p^.^+i ) ( n ^0 ( n h.^+i ] (^ = -i ^my) . (12) 

Although the construction f lT2|) takes advantage of the sparse matrix decomposition, it 
is not maximally efficient for reasons that we now point out. 

Consider first the case of generic v. Successive applications of T to |eo), as in 
(E]), will generate a number of basis states corresponding to all possible partitions of 
{1,2, . . . ,m} that are consistent with planarity. These are known in the literature as 
non-crossing partitions, and their number is given by the Catalan number 

^- = ,P"l'n, = 1' 2, 5, 14, 42, 132, . . . , (13) 

m\[m + Ij! 

which grows asymptotically as Cm ~ ■m~'^/^4™' for m 3> 1. 
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For V = —1 we can use (IT^ . and so once a layer is completed the number 
of basis states equals the number of non-crossing non- nearest- neighbour partitions of 
{1,2, ... ,m}. Their number is given [3] by the Motzkin number Mm-i, where 

^^= E (2-)c, = l,l,2,4,9,21,.... (14) 

For instance, when m = A, the four possible such partitions read 

(1)(2)(3)(4), (1)(24)(3), (13)(2)(4), (14)(2)(3). 

The asymptotic behaviour is now M^ ~ /c^^/^3'^ for A; ^ 1. This looks like a considerable 
improvement. However, the intermediate stages that appear when building up the 
product flT2|) still require all nearest-neighbour partitions, and so the time necessary for 
a complete multiplication with T still grows like 4™. 

Fortunately, we can write the sparse matrix decomposition of T in a way which 
(essentially) only makes use of non-crossing non-nearest-neighbour partitions. This is 
accomplished by utilising the property ([9]) maximally. Indeed we can write, for v = —1, 

T = T,nTm-i---nT2T, (15) 



with 



m—l 



ViHi,i+i for i = 1,2,. . .,m - 1 






(m—l \ 

llP,,J+lj'^m. (16) 

After the multiplication by Tj, the non-nearest-neighbour constraint is now in action 
everywhere, except possibly between sites i and i + 1. The number of partitions where i 
and i + 1 do belong to the same block is Mm-2, since for such partitions the contraction 
of i and i + 1 into a single point clearly leads to a non-crossing non-nearest-neighbour 
partition on ?7i — 1 points. 

The sparse-matrix factorisation (TT5l) - (TT6l) is a central result of this paper. It implies 
that, for the square lattice, the number of basis states needed to compute the chromatic 
polynomial using IQ and flTSl) is M^-i when a layer has just been completed, and 
Mm-i + Mm-2 in the intermediate stages. Both numbers grow like 3™. 

It remains to describe how the matrix Tiast lays down the last layer of horizontal 
edges. As before this is done one edge at a time 

last - J-m ^m-1 ' ' '^3 ^2 ^1 U ' J 

In this product, each factor adds a horizontal bond between i and i + 1. Having done 
that, vertex i and those to the left of it are of no use any longer, so we might as well 
detach them as singletons. This leads to 

^last ^ I JJ p^.^.^^ j [)^H,^,+i for z = 1, 2, . . . , m . (18) 
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Note that after the action of Tiast, only the all-singleton state |eo) remains. The 
chromatic polynomial of G appears as the coefficient of this state, leading finally to 

3.3. Triangular lattice 

The triangular lattice can be drawn as a square lattice with diagonals. These diagonals 
can be constructed, within the transfer matrix approach, by applying the operators 
Hj_j+i just after site i has been propagated from one layer to the next by means of Vj, 
but before Vj+i has been employed. 

We have again the factorisation ( ITSjl . but now with 

(m-l \ 

/m-l \ 

T;; = I TT P,.i+i I yMi-i,Mi,i+i for i = 2, 3, . . . , m - 1 

Tm = I J_J_ PjJ+l I ^m^m-l,m- (19) 

Note that, due to the existence of the diagonal edges, the leftmost factor in each 
expression is the product over all nearest-neighbour projectors. Therefore, the number 
of basis states required in any intermediate stage of the construction is Mm-i- 

The expressions ( 1T7|) -( 1T8|) for Tiast are identical to those for the square lattice. 

3.4. Ordering the states 

At this point, one would be ready to write a transfer matrix for the chromatic 
polynomial using only non-crossing non-nearest-neighbour (NCNNN) partitions — up to 
a small subtlety for the square lattice — as intermediate states, by storing the partitions 
encountered using standard hashing techniques. 

One can however reduce further — albeit only by a constant factor — time and 
memory requirements by establishing a bijection between the natural numbers 
{1,2, . . . ,Mm-i} and the set of M^-i NCNNN partitions. Using this, the vectors 
appearing in the computation of (jH]) can be realised simply as standard arrays, indexed 
by the natural numbers, and each basis state can be mapped back and forth between its 
representation as a natural number and as a partition on which the operators Hj_j+i, Vj 
and Pi,i+i can be made to act. This idea was previously implemented in the context of 
the planar Potts model (i.e., for non-crossing partitions) by Blote and Nightingale [16], 
and we extend it here by adding the non- nearest-neighbour constraint. 

An important step in establishing the bijection is to produce an ordering of the 
NCNNN partitions for a given m. To this end, consider the first point in the partition. 
Two situations may arise: 
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(i) If point 1 is a singleton, the restriction of the partition to the points {2, 3, ... , m} 
is again NCNNN. 

(ii) If point 1 is not a singleton, let A; G {0, 1, . . . , m — 3} be the smallest integer so that 
point 3 + k belongs to the same block as 1. Then the restriction of the partition to 
points {2,3, ... ,2 + k} and the restriction to points {3 + fc, 4 + /c, . . . , m} are two 
independent NCNNN partitions of respectively k + 1 and m — k — 2 points. 

This argument proves in particular the recursion formula (with p = m — 1) 

p-2 

Mp = Mp_i + J2 MkMp_k.2 (20) 

fc=0 

with initial values Mq = Mi = 1. 

Suppose now, as an induction hypotheses, that we have defined an ordering of 
NCNNN partitions of fewer than m points. The ordering of NCNNN partitions of 1 or 
2 points is trivial, since in these cases there is only one such partition. We can then 
order NCNNN partitions of m points as follows: 

(i) Place first the partitions in which point 1 is a singleton. The ordering of 
these partitions follows inductively from the ordering of the partitions of points 
{2,3,...,m}. 

(ii) Place next the partitions in which point 1 is not a singleton, in the order of 
increasing values of the integer k defined above. The ordering of the partitions 
with a fixed value of k follows inductively from the ordering of the partitions of 
points {2, 3, . . . , 2 + A;} and {3 + k,4: + k, . . . , m}, using the former (resp. latter) as 
the most (resp. least) significant 'bit'. 

For instance, the ordering for ?7i = 5 points reads: 
(1)(2)(3)(4)(5) , (1)(2)(35)(4) , (1)(24)(3)(5) , (1)(25)(3)(4) , 
(13)(2)(4)(5) , (135)(2)(4) , (14)(2)(3)(5) , (15)(2)(3)(4) , (15)(24)(3) . 
The mapping from integers to partitions can now be constructed by inferring the 
value of k from (l20l) and proceeding recursively. The inverse mapping is constructed 
similarly by obtaining k explicitly from the partition. 

Finally, for the square lattice, we need a bijection between the integers 
{1,2, . . . , Mm-i + Mm-2} and partitions in which points i and i + 1 are allowed to 
be in the same block. Again this follows from our ability to order such partitions. To 
this end, we place first the M^-i partitions in which i and i + 1 are not in the same 
block, and next the Mm-2 partitions in which i and i + 1 are in the same block. The 
remaining ingredients of the bijection can be taken over directly from the NCNNN case. 

3.5. Implementational details 

We wish to compute the chromatic polynomial, normalised as in ([7]), 

PG(g) = g~i''i^G(g,-i), (21) 
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where G is an m x n rectangular section of the square or triangular lattice (the latter 
being considered as a square lattice with added diagonals). The result is a polynomial 
of degree \V\ — 1 in the variable q~^ with constant term 1. To obtain this by the 
transfer matrix method, each NCNNN state must have a coefficient (weight) which is 
a polynomial of degree d = |\^| — 1 with integer coefficients. These coefficients can be 
stored quite simply as an integer array of length d+ 1. 

Since the goal is to produce a series expansion in powers of q^^ of the various free 
energies, there is actually no need to compute the whole polynomial Pcil)- We shall 
see below that one knows in advance to which degree (imax the expansion can be taken. 
It suffices therefore to truncate all polynomials intervening in the coefficients to that 
order. Such a truncation leads to non-negligible savings of memory resources. 

A practical problem arises because the integer coefficients rapidly become very 
large. This is effectively handled by using modular arithmetic, i.e., by computing 
the coefficients Pg{(i) modulo various primes Pi. Since the transfer process only uses 
additions and subtractions, it is convenient to use 32-bit unsigned integers in the 
computations and to choose Pi < 2^^ in order to avoid overflow before each modulo 
operation can be carried out. Using the Chinese remainder theorem, the partial results 
Pg{q) mod Pi can be combined to give Pail) mod HiPi- Once the latter stabilises upon 
using further primes we know that we have arrived at the complete result Pciq), with 
no modulo. The results reported below made necessary the use of up to 16 different 
primes. 

4. Results 

We first review briefly how the series expansion for the free energies can be obtained 
from the chromatic polynomials on flnite m x n rectangles Gm,n by use of the flnite- 
lattice method [17]. The goal is to obtain a development for the free energy per vertex 
in the form 



Jm.,n\C[) 



logPciq) 



mn 



(22) 



log g + ^ 



fe=0 



Ofc 

gk 



+ 



.fc=0 



1 1 
— + - 

m n 



+ 



ECk 
qk 

.k=0 ^ 



1 

mn 



exactly valid, order by order in g ^, for all m and n sufficiently large. The quantities 
within square brackets are referred to as the bulk, surface and corner free energies, fh{q), 
/s(g) and /c(g). 



4.I. Finite-lattice method 

We discuss first the case of the square lattice. The free energies introduced above are 
given by the flnite-lattice method [17] as 

/b(g) = logg + J2 "'^('^' ^) ^ogPG(r,s){^) + Oiz"'-') , (23) 

{r,s)eB{k) 
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fs{q)= Yl Pk{r,s)\ogpGir,s){z) + 0{z''-'), (24) 

{r,s)eBik) 

Uq)= Y. lk{r,s)\ogPGir,s){^) + 0{z^^-^), (25) 

(r,s)£B(k) 

where the set B{k) are r x s rectangles G{r,s) such that r < s and r + s < k. 
For convenience — and to facihtate the comparison with previous work [181 IH] — the 
expansion variable q~^ has here been changed to 

z = ^. (26) 

q-l 

The error term 0{z'^^~^) in (123|) - (!25|) is due to certain classes of convex polygons [T7] 
which do not fit into any of the rectangles in B{k). 
The finite-lattice weights for /b are given by [IT] 

Q^fc(^, ■S) = (2 — 6r,s) {Sr+s,k — 35r+s,fc-l + 35r+s,A:-2 " (^r+s.fc-s) (27) 

while those for /g and /c read [21j 

fl r ^ ^9 A ^ /^2-r-s 3(r + s)-2 

Pfc(r, s) = (2 - ^r,J I dr+s,k H ^ d^,+s,fc_i- 

3(r + s) + 2 , r + s-2 ^ \ 



2 ■-=-- '2 

7A:(r, s) = (2 - 5r,s)((r - l)(s - l)(5r+s,fc + (1 + r + s - 3rs)6r+s,k-i + 

(3rs + r + s - l)(5,+,,fc_2 - (r + 1) (s + l)6r+s,k-3) ■ (29) 

-^.^. Square lattice series 

Using the transfer matrix formalism we can compute PG{r,s) for lattice strips of width r < 
^max- Without resorting to super computer facilities or massive parallel computations, 
we have found that mmax = 20 was feasible on a standard workstation. This corresponds 
to a maximum dimension of the transfer matrices of Mig + Mig = 24 735 666 on the 
square lattice. Setting fc = 41 above, this gives then the free energies /b, /s and /c 
correctly to order z'^^. For the case of fy, and /g an extra term can be obtained by adding 
manually the term fi2mm.^z^"''"""~^ to /b [E], and the term -mmaxyU2m,„ax^^"'"""~^ to /s 
[11], where [M] 



V- 2p_ 8 2-2x^-xVl-4x2 i2 3-4x^-4v/l-4x2 

^^^"^ "^ (1-4X2)(2 + X2) +^ (1-4x2)2 • ^"^"^ 

We thus obtain the series to order z^^ in these cases. 

The results for all three series are given in Table [H They extend those of [H] by 
32 new terms. 

It is possible to analysq^l the extended series using the method of differential 
approximants. The goal is to obtain the locations Z£ and exponents A^ of the singularities 

I We are grateful to J. Salas for providing the analysis mentioned in this paragraph. 
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14 


8 


-210 
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47 
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849 
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1471 
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254 


-493 


10020 
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1671 


-8052 


30548 
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3127 


-19901 


68832 
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786 


-19966 
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-13939 


37556 


65229 
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-49052 


223807 


-236055 
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-80276 


508523 


-739289 
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21450 


321314 


101404 


26 


515846 


-2052462 


7201383 


27 


1411017 


-8417723 


26255714 
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1160761 


-13374892 


43505098 


29 


-4793764 


10841423 


-17552274 


30 


-20340586 


112595914 


-291420026 


31 


-29699360 


260687001 


-674637832 
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33165914 


70989018 


27442 


33 


256169433 


-1341964856 


4426763291 


34 


495347942 


-4108283969 


12910062402 



Table 1. Large-g series for the bulk, surface and corner free energies of the chromatic 
polynomial on the square lattice, in terms of the variable z = l/(<7 — 1). We 



give the coefficients of exp(/b) 



_ jq-i) 



J2kLo(^kz'', exp(/s) = Efclo/^fc^'"' and 
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k 


ttfc 


^k 


Ik 


35 


-127736296 


-3304416038 


9737827437 


36 


-3068121066 


14960606999 


-49131891078 


37 


-7092358808 


58237169596 


-184470253912 


38 


-1024264966 


65268280922 


-183956055539 


39 


35697720501 


-162368154719 


621518352427 


40 


91243390558 


-767619757924 


2660084937207 


41 


25789733672 


-975329692910 


3075466954690 


42 


-420665229170 


1872486336165 


-7500763944932 


43 


-1089052872105 


9701425034093 


-34822638005931 


44 


-238516756366 


12262136381593 


-38841312202313 


45 


5101697398582 


-24192583755347 


104412348649015 


46 


12146149238921 


-118764516172484 


448320847685638 


47 


-598931311074 


-130312353695974 


421063156900936 


48 


-63338329084478 


346985943489639 


-1507607631713074 


49 


-125863684143541 


1412474002216034 


-5446631048235906 


50 


71523371777335 


1098760305055003 


-3234717862390686 


51 


790320197578279 


-5183269098276812 


22683201852542681 


52 


1200257719380100 


-16304598291388925 


64013086293590143 


53 


-1721302814683702 


-5448655170720498 


8336294391279875 


54 


-9694394826555237 


75813721738614619 


-329171853716322290 


55 


-10463394611604378 


182868102163636571 


-718400221457685146 


56 


29575960135787439 


-33129258404821172 


310674267788777570 


57 


115472492262862427 


-1051108217408277081 


4563895897602278816 


58 


85395764100885186 


-2017022924530473193 


7953356311623456864 


59 


-421383549369627730 


1373959176501915516 


-7763050320815504003 


60 


-1343046574902992458 


13749175266084267367 


-60017851952733958335 


61 


-736061882807769534 


22595026341744522020 


-90527220474610273313 


62 


5291983916641725051 


-23294356475365888405 


121713573010967365453 


63 


15613586538385554242 


-172293246256852139137 


768987768617382802911 


64 


8244136015720127003 


-268117182162694379708 


1123758858144026736112 


65 


-61540273845219665190 


298116598787375322521 


-1553871017012870302628 


66 


-186576584391637858059 


2127014580990734112356 


-9884276597620937431531 


67 


-119156912656052312977 


3423100771562123299025 


-15168542715543012752367 


68 


700439920487427255618 


-3325599173839179924758 


18496951095018457216592 


69 


2311567965527717105415 


-26543042275747617969190 


130030721001519590513331 


70 


1784944419294044372980 


-45519479688835787926364 


210058368547617077899237 



Table 1. (continued. 
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k 


ftfc 


/3fc 


Ik 


71 


-8190295787656754694494 


36371250229750107503644 


-228752490963172501272677 


72 


-29120128759101310463081 


336827820960866947369728 


-1729404450220843332781816 


73 


-24509224734165974718725 


600304759572952426851 130 


-2809877317417695436535077 


74 


99609960301714320809469 


-429784892999554830063487 


3052097912541295946004847 


75 


363523561198034374612305 


-4289105251916055148255439 


22691488586148969734635242 


76 


307076608628563447831354 


-7658381117547376351931511 


35816140480454963887649892 


77 


-1234243320099804373986722 


5508060829297214833104965 


-41473851023699412644120745 


78 


-4471768398561919868797416 


54201975484843398706779272 


-291274582715872945456529116 


79 


-3688741741221487636404873 


95450369797304235548078036 





Table 1. (continued. 



of the free-energy series, which is supposed to behave hke ~ (z — -2^)^* as z — )■ zi. 
Proceeding exactly as described in section 5.5.3 of [H], one obtains for the bulk free 
energy the first-order {K = 1) approximants 

zi = 0.24(4) ± 0.40(5)i Ai = -1.8(7) 

Z2 = 0.64(7) As = -1.6(7) (31) 

23 = -1.00(1) A3 = 0.99(6) 

while the second-order {K = 2) approximants read 

z, = 0.24(4) ± 0.40(3)2 Ai = -2.4(7) 

Z2 = 0.62(10) A2 = -1.3(16) ^ ^ 

Somewhat disappointingly, these results hardly improve on those reported in |llj . 
despite of the substantial extensions of the series. The same is true for the reanalysis of 
the surface and corner free energies. 



4-3. Triangular lattice series 

For the triangular lattice, we can use exactly the same finite-lattice weights as for 
the square lattice. Taking again mmax = 20 and k = 41 corresponds to a maximum 
dimension of the transfer matrices of M19 = 18 199 284. The error term in ( |23l) - (!25ll is 
now 0(2^^), and we know of no special tricks for obtaining an extra term. Therefore, we 
obtain in this case the three series to order z^^ included. 

In a couple of important papers, Baxter [191 E] has shown that the triangular- 
lattice chromatic polynomial is an integrable system. In particular, he has set up the 
Bethe Ansatz equations and derived three different analytical expressions for the bulk 
free energy /b. The question which one of the three expressions dominates for a given 
real — or even complex — value of q turns out to be very intricate [6], and has recently 
been the subject of a careful reanalysis p]. The expression which is dominant for real 
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q G (—00, 0] U [4, 00) turns out to be 

X y (1 - a;6j-5)(i _ a;6j-4)(i _ a;6j)(i _ ^6j+i^ ^ > 

where the previous expansion variable fl26l) has now been traded for x which is defined 
by 

g = 2 - X - - . (34) 

X 

As expected, our series expansion vahdates (!33|) to order x^^ included. 

The analytic computation of /s and /c was not attempted in [T^l E]. Presumably 
/s could be derived by a substantial adaptation of Baxter's work in order to impose 
non-periodic boundary conditions in his transfer matrix formalism. On the other hand, 
the derivation of /c appears to be a much harder problem, since this would require 
precise knowledge of the correspondence between boundary states in the Bethe setup 
and the initial lattice model, as well as the ability to compute scalar products in the 
Bethe formalism. 

The series expansions that we obtain for /s and /c read as follows 

e^= = 1 - 2x + x^ + 2x^ - 2x^ - x^ + 2x^ + 2x^ - 4x^ - 2x" + 7x^^ - 
8x^^ + 4x^^ - hx^"^ + 12x^^ - 12x^^ + 4x^° - lOx^^ + 24x^^ - 
2x^^ - 19x^^ + 2x^^ - ISx^^ + 46x2^ - lOx^^ - 24x^^ - x^° - 
24x^^ + 77x^2 - 22x^^ - 25x^^ - 16x^^ - 35x^^ + 118x^^ - 
41x^^ - 14x^9 - 45x^° + 0(x^^) , (35) 

e^- = 1 + X - 3x^ - 5x^ - 6x^ + 2x^ + llx^ + 21x® + Ox^ - 12x^° - 

41x^^ - 33x^2 + 72x^^ + 82x^^ + 40x^'^ - 102x^^ - 169x^^ - 

139x^^ + 103x2° + 273x2^ + 319x^2 - SOx^^ - 398x2^ - 612x2^ - 

lllx^*^ + 517x2^ + 1066x2^ + 454x2^ - 537x^° - 1686x^^ - 

1067x^2 + 398x^^ + 2500x^^ + 2079x^^ + 76x^'^ - 352 Ix^^ - 

3641x^^ - 1138x^9 + 4643x^° + 0(x^^) . (36) 

It is seen that the coefficients are numerically much smaller than those for the square 
lattice (see Tabled]). This in itself is not surprising, since the radius of convergence is 
1 in the x variable (triangular-lattice case), while that of the z variable (square-lattice 
case) is ~ I (and presumably exactly | in view of [26j). What is truly surprising is 
that the coefficients are integers! It thus looks reasonable to try to conjecture exact 
expressions. 

In analogy with other problems where both bulk and surface properties have been 
worked out, one would expect the periodicity of 6 in the product formula ( l33ll to be 
doubled to 12 in the corresponding surface quantities. Taking this as a clue suggests 

trying for e-^'' and e-^'^ product formulae of the form 

00 11 
JJJJ(l-x^'^-'=)°'= (37) 

j=\ fc=o 
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, 12 

5 ■' 



for some exponents a^ to be determined. Knowledge of the coefficients oi x^, . . . ,x 
from fl35|) -f l36|) allows us to determine ao,...,aii; we can then test our conjectural 
formula against the coefficients of x^^, . . . , x^^. 
In this way we have obtained the conjecture 

y V (1 - a;i2j-9)(l - xi2j-8)(i _ a;i2i-7)(i _ 3,12^-2) J ^ ) 

for the surface free energy, and 

00 

e^^ = 11(1 - x''^-'Y\l - a;i2^-i°)(l - x'^^-y{l - x'^^-^'f x 
j=i 

(1 _ a;i2^-7)4(l - a;i2^-6)-3(i _ xi2^-5)(l - a;^2i-4)-3 ^ 

(1 _ a;i2^-3)5(i _ a;i2,-2)2(^ _ xi2^-i)2(l - x''^ (39) 

for the corner free energy. The validation of fl551) - fl5I?l) by 40 — 12 = 28 independent 
checks means that both conjectures are established beyond any reasonable doubt. 

4.4- Limits g — )■ and g — )■ 4 

It is obvious that the product formulae fl33|) . fl38|) and fl39l) converge for |x| < 1, since 
this is true for each individual factor in the product (by taking logarithms). It remains 
to discuss the limiting cases x — )■ 1 (i.e., g — )■ 0) and x — )■ — 1 (i.e., g — )■ 4). The 
latter case — four-colouring the vertices of the triangular lattice — can be shown to be 
equivalent to other physically interesting problems, such as three-colouring the edges 
of the hexagonal lattice [221 123], or fully-packed loops on the hexagonal lattice with 
fugacity nioop = 2 [24]. 

In the X — )• 1 case, we write x = 1 — e for e ^ 1. Then, to leading order, we have 
for each factor in the products 1 — x" ~ ae. This leads to 

(6j-3)(6j-2)2(6j-l) 



lim e-^'' = — TT 

rr^^ J- J- 



i-z;i ' LL ^Qj _ 5) (6j - 4) (6j) (6j + 1 

33/^r (If 



(27r)5 



3.770919 ■ ■ ■ , (40) 



X— >1 



lime-'^ = I I 

J = l 
1 

~ 2^ 



• (12j-ll)(12j-6)(12j-5)(12j-4) 
. (12j-9)(12j-8)(12j-7)(12j-2) 

r(^)rQ) 



2 



_0.305639--- , (41) 

r(|) 

hm e-''^ = . (42) 

In the X — !■ — 1 case, we write x = —1 + e for e ^ 1. When a is even, we have then 
x" ~ ae as before. And when a is odd, 1 — x"" ~ 2 to leading order. This leads to 

lim e^'' = n ,}.^ ~ L., = }l' ^ 1.460998 ■ ■ ■ , (43) 

"^-^ |i(3j-2)(3j) 47r2 ^ ^ 
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lim e-''' = I I 



i=i 



f'~T''';\ ] ' = 4^/^ - 1-587401 • • • , (44) 

_(6j -4)(6j - 1)J 

lim e-''' = oo . (45) 

It is remarkable that /c diverges in both limits x — )■ ±1. 

5. Discussion 

In this paper we have used the finite-lattice method to obtain large-g series for the bulk, 
surface and corner free energies — denoted /b, /s and /c — for chromatic polynomials 
on regular two-dimensional lattices. Long series could be obtained due to a technical 
improvement of the transfer matrix method, which is specific to the antiferromagnetic 
zero-temperature limit of the Potts model, and amounts to obtaining a factorisation of 
T in terms of sparse matrices, each of dimension ~ 3"^, where m is the strip width. 
Although this applies in principle to any regular lattice, we have focussed here on the 
most commonly studied lattices (square and triangular). 

For the square lattice the improvement allowed us to add 32 new terms to the 
existing series [11] . Unfortunately, the reanalysis of the singularities of this series, using 
the method of differential approximants, did not lead to any significant improvement 
over [H]. The focus in this paper was therefore rather on the triangular lattice, where 
the series permitted us to conjecture the exact analytical formulae (138|) -( 139|) for /g and 
/c, and to confirm the correctness of Baxter's [19], E] result ( l33l) for /b. 

It should be noticed that going to mmax = 20 leads to 79 terms in the square-lattice 
series, but only to 40 terms for the triangular lattice. Longer series could be obtained in 
the latter case by using a version of the finite-lattice method which is specifically adapted 
to the triangular lattice [25]. However, since our series have enabled us to conjecture 
the exact (infinite-order) results (!38l) -(l39l ) the motivation for extracting further terms of 
the series is admittedly rather scarce. Indeed, the series reported here already provide 
40 — 12 = 28 independent checks of the proposed expressions ( l38l) -(!39l). 

It would be interesting to elucidate the physical mechanism by which f^ diverges 
in the limits g — )■ and g — )■ 4. 

It is regrettable that the series obtained here do not provide any information about 
the regions where the chromatic polynomial is a critical system — that is, q G [0, 3] for 
the square lattice [26l[l3l[Tll|3l[7l[in] and g G [0, 4] for the triangular lattice [T9| |6| |8| [TO] . 
The reason is that g = oo seems to be the only viable expansion point for the finite- 
lattice method. 

It is possible that the fact that non-nearest neighbour constraint can be made 
local conceals some deeper representation theoretical meaning. It is well-known that 
the generic representation theory of the Temperley-Lieb algebra, which is central to 
the algebraic formulation of the planar Potts model, undergoes profound modifications 
when g = 4cos(7r/(p + 1))^ with p integer (or rational) [27J. It might be that similar 
modifications occur upon taking the particular value v = —1 for the temperature 
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variable. It should be noticed that some steps toward an algebraic formulation of the 
chromatic polynomial have been taken in 
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